Astronomy & Astrophysics manuscript no. 5220 


© ESO 2008 


February 5, 2008 





o 
o 



> 



The Gamma-ray burst 050904 : evidence for a termination shock ? 

B. Gendre 1 , A. Galli 1 - 2 ' 3 , A. Corsi 1 ' 2 - 4 , A. Klotz 5 - 6 , L. Piro 1 , G. Stratta 7 , M. Boer 6 , and Y. Damerdji 6 - 5 



IASF-Roma/INAF, via fosso del cavaliere 100, 00133 Roma, Italy 

e-mail: bruce . gendre@iasf-roma . inaf. it , alessandra . corsi@iasf-roma . inaf . it , galli@iasf-roma . inaf . it , 

piro@iasf-roma . inaf. it 

Universita degli Studi di Roma "La Sapienza", Piazzale A. Moro 5, 00185, Roma, Italy 
, 3 INFN - Sezione di Trieste c/o Dipartimento di Fisica - Universita di Trieste Via Valerio 2, 34127 Trieste, Italy 
£ — ' 4 INFN - Sezione di Roma c/o Dipartimento di Fisica - Universita degli Studi di Roma "La Sapienza" Piazzale A. Moro 5, 00185 Roma, Italy 
(N ', 5 CESR, 9 avenue du Colonel Roche, 31400 Toulouse, France 
e-mail: klotz@cesr.fr 

Observatoire de Haute-Provence, 04870 St Michel l'Observatoire, France 
e-mail: michel . boer@oamp . fr 

LATT, Observatoire Midi-Pyrenees, 14 avenue Edouard Belin, 31400 Toulouse, France 
e-mail: gstratta@ast.obs-mip.fr 



m 
m 

Received — ; accepted — 

\o : 

O ■ ABSTRACT 

O | Aims. We analyze optical and X-ray observations of GRB 050904 obtained with TAROT and SWIFT. 
I ■ Methods. We perform temporal and spectral analysis of the X-ray and optical data. 

O Results. We find significant absorption in the early phase of the X-ray light curve, with some evidence (3<x level) of variability. We interpret this 
■4_> as a progressive photo-ionization. We investigate the environment of the burst and constrain its density profile. We find that the overall behavior 
^ ' of the afterglow is compatible with a fireball expanding in a wind environment during the first 2000 seconds after the burst (observer frame). 
. . On the other hand, the late (after 0.5 days, observer frame) afterglow is consistent with an interstellar medium, suggesting the possible presence 
^ of a termination shock. We estimate the termination shock position to be R, ~ 1.8 x 10~ 2 pc, and the wind density parameter to be A, ~ 1.8. 
' We try to explain the simultaneous flares observed in optical and X-ray bands in light of different models : delayed external shock from a 

thick shell, inverse Compton emission from reverse shock, inverse Compton emission from late internal shocks or a very long internal shock 
g$ activity. Among these models, those based on a single emission mechanism, are unable to account for the broad-band observations. Models 

invoking late internal shocks, with the inclusion of IC emission, or a properly tuned very long internal shock activity, offer possible explanations. 
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1. Introduction 

Long Gamma-Ray Bursts (GRBs) are cosmological explosions 
(Metzger et al. 1997), thought to be produced by a massive 
star (MacFadyen & Woosley 1999). Their extreme luminosity 
(up to L, J0 ~ 10 53 erg s _1 , Piran 2005) makes them observable 
up to very large distance (z> 5). Their long lasting afterglow 
allows one to study them for days after the trigger (for a recent 
review see Meszaros 2006). 

With the launch of the SWIFT satellite (Gehrels et al. 

2005) , it has been shown that strong flaring activity in the early 
X-ray afterglow originally found in a few cases by BeppoSAX 
(e.g. Piro et al. 2005) is a common phenomenon (O'Brien et al. 

2006) . The origin of these X-ray flares is not clear, and several 
models have been proposed: late internal shocks (Fan & Wei 



2005; Burrows et al. 2005b), refreshed energy injection due to 
long lasting activity of the GRB progenitor (Zhang et al. 2005), 
rising of an extra spectral component such as inverse Compton 
(Kobayashi et al. 2005), late internal and external shocks (Wu 
et al. 2006), delayed external shock emission in a thick shell 
fireball (Piro et al. 2005; Galli & Piro 2006). While some 
of these flares show spectral evolution and may be linked to 
the prompt emission (Burrows et al. 2005b), others look very 
similar to the late afterglow emission and can be linked to its 
onset (Piro et al. 2005; Galli & Piro 2006). Nevertheless, no 
multi-wavelength observations are usually available to discrim- 
inate these hypothesis, and several models can fit the data (e.g 
GRB 050406, Wu et al. 2006). The time dilation, due to cos- 
mological effects, mostly evident for high-z bursts, allows to 
study more easily the early afterglow and the prompt-to after- 
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glow transition. Bursts located at high distances are thus very 
interesting to address these issues. 

Another interesting aspect of the afterglow studies is the 
constraints one can put on the absorption around the burst 
(Stratta et al. 2004; Campana et al. 2006). The absorption 
in the X-ray band is primarily determined by metals. Using the 
low energy part of the X-ray spectrum, one can thus constrain 
the metal content of the burst surrounding medium. This is a 
key issue for high-z bursts, as it allows to study the metal en- 
richment of the Universe at early epochs. 

A final aspect of the afterglow evolution is the problem of 
the surrounding medium of the burst. Since long GRBs are re- 
lated with massive stars, one should expect a wind environment 
around them (Chevalier et al. 2004). On the other hand, several 
GRBs exhibit afterglows consistent with an expansion in a con- 
stant density medium (hereafter InterStellar Medium, or ISM) 
rather than a wind environment (Panaitescu & Kumar 2001). 
Several authors (e.g. Ramirez-Ruiz et al. 2001) proposed that 
the expanding wind arising from the star can be stopped by a 
dense surrounding interstellar medium. The interface between 
these two media is called the termination shock. In this sce- 
nario, we would expect the fireball expansion to be consistent 
with a wind only in the early phase of the afterglow, and con- 
sistent with an ISM later (Chevalier et al. 2004; Ramirez-Ruiz 
et al. 2001). However, such a transition was never observed so 
far. Here we present a multi-wavelength analysis of the GRB 
050904 afterglow and show its consistency with the occurrence 
of a termination shock. 

GRB 050904 was at z= 6.29 and presented strong and long 
lasting multi-wavelength flaring activities (Kawai et al. 2006; 
Watson et al. 2006; Boer et al. 2006). We summarize previous 
observations in Sec. 2 and describe the data analysis in Sec. 
3. We report in particular a more detailed analysis of the X- 
ray and TAROT data. We investigate the absorption evolution 
around this burst in Sec. 4, and the nature of the progenitor to- 
gether with the medium density profile in Sec. 5. We discuss 
the origin of the observed multi-wavelength flare in the frame- 
work of a delayed external shock emission due to a thick shell 
fireball in Sec. 6.1. Other models proposed in the literature are 
discussed in Sec. 6.2. Conclusions are given in Sec. 7. 

2. GRB 050904 

GRB 050904 was triggered by the BAT instrument (Barthelmy 
et al. 2005) on board the SWIFT satellite (Gehrels et al. 2005) 
on September 4th, 2005, at 01:51:44 UT (Cummings et al. 

2005) . This was a bright burst, with a fluence of 5.4 + 0.2 x 10~ 6 
erg crrT 2 (15-150 keV), and a duration of Tgo = 225 + 10 sec- 
onds (Sakamoto et al. 2005). The energy spectral index was 
0.34 + 0.06 (Sakamoto et al. 2005). The narrow field instru- 
ments XRT (Burrows et al. 2005a) and UVOT (Roming et al. 

2006) observed the field of GRB 050904 about 160 seconds af- 
ter the trigger. No optical counterpart was observed within the 
UVOT field of view, while an X-ray transient was detected at 
position 00* 54 m 50.4 s +14° 05' 08.5" (Cummings et al. 2005). 
The first optical detection was made by TAROT (Boer et al. 
1999), which observed a faint optical afterglow in the unfiltered 
frames (Klotz et al. 2005b,c; Boer et al. 2006). Other observa- 



tions made with the BOOTES-IB telescope failed to detect this 
afterglow in the R band (Jelinek et al. 2005). Larger optical 
and infrared telescopes imaged the field of view 3 hours after 
the burst and detected an infrared counterpart (Haislip et al. 
2006). The absence of any counterpart observed from the R to 
U bands implied a large absorption, possibly due to the Lyman 
alpha forest (Haislip et al. 2006). This was confirmed by the 
spectroscopic redshift measured by the Subaru telescope (z = 
6.29 Kawai et al. 2006), which implied that the Lyman alpha 
cut-off was redshifted in the optical up to 8800 A. This burst is 
currently the more distant burst ever observed. 

3. Data reduction and analysis 

3.1. X-ray data 

3.1 .1 . Data reduction 

We obtained the XRT data from the SWIFT archive 1 and re- 
duced them using the available packages and calibration files 
(SWIFT FTOOLS version 2.2 and CALDB version 2005 1028). 
Data were filtered using the provided good time intervals (sun 
and moon constraints, bright earth and South Atlantic Anomaly 
limitations) and standard criteria : CCD temperature < -50°C; 
grades to 2 and energy above 0.3 keV for Window Timing 
mode (WT); grades to 12 and energy above 0.5 keV for 
Photon Counting mode (PC). Note that the spectral ranges are 
different in WT and PC mode in order to take into account some 
uncertainties on the PC mode calibration (Osborne et al. 2005). 
We extracted spectra and light curves using circle (PC mode) 
or box (WT mode) regions of 25 pixels radius (that enclose 
~ 95 % of the PSF, as explained in the XRT data reduction 
guide 2 ). The background was estimated using a larger region 
free of sources. 

Our filtering criteria are mode dependent due to the change 
in the spectral range, and the mode switch may introduce 
some problems due to cross calibration uncertainties. Thus, any 
count light curve extracted need to be corrected for this depen- 
dency. To do so, we constructed a color-color diagram to look 
for spectral variability. This information allowed us to convert 
the count light curve into a flux one which is mode indepen- 
dent. This light curve is presented in Fig. 1. As one can see, 
the transition between the two modes occurs at the end of a 
large flare; however, the continuity of the light curve made us 
confident on our analysis. 

We extracted spectra for several time intervals listed in 
Table 1 . The spectra were re-binned to contain at least 20 net 
counts per bin using GRPPHA, and were fitted using XSPEC 
version 1 1.3.1 (Arnaud 1996). All errors for the fit parameters 
quoted in the paper are given at the 90 % confidence level for 
one interesting parameter. 

3.1 .2. Data analysis 

The X-ray light curve shows a very complicated evolution, with 
several flares. In the following, we take as the start of a flare 

1 see http://swift.gsfc.nasa.gov/docs/swift/archive/ 

2 available online at http://swift.gsfc.nasa.gov/docs/swift/analysis/ 
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the point from which the derivative of the decay law becomes 
positive. The end of a flare is defined as the point where the 
light curve shows a flattening (we consider the flattening as the 
recovery of the continuum level). The early X-ray light curve 
(t < 2000 s) shows a steep decrease, followed by a plateau 
and a large flare (see Fig. 1). A second flare occurs at 1240 
seconds (Mineo et al. 2005). We fitted this part of the light 
curve excluding all data points between the start and the end of 
all flares. Since the estimation of the start date of the first flare is 
complicated by the presence of the plateau, we conservatively 
excluded also the data located in the plateau (see Fig. 1). Since 
a sudden softening of the spectrum is observed at the time of 
the first flare, we have fit separately the data before and after 
the flare using a power law. The early decay index is 2.5 ± 0.2 
ixl = 1.05, 13 d.o.f.) and the late one 1.5 + 0.2 (x 2 Y = 0.72, 14 
d.o.f.). 

The whole spectra (0.3 keV or 0.5 keV to 10 keV) were fit- 
ted with a simple power law (i.e. F(v) oc v~ a ), absorbed by our 
galaxy (Nh value fixed to the galactic one of 4.97 x 10 20 cm -2 
Dickey & Lockman 1990) and by local absorbers at redshift 
6.3. We detect an excess of absorption before 264 s. After that, 
we have only upper limits, except between 384 and 463 sec- 
onds (which corresponds to the flare maximum). In some cases, 
the XRT may experience some gain variation for bright sources 
due to a bias subtraction problem (Campana et al. 2006). In 
such a case systematic deviations from the model appear near 
the oxygen edge (0.5-0.6 keV). Excluding the spectral bins be- 
tween 0.5 and 0.6 keV from the fit of the flare spectrum, the 
Nh is found to be 2.4^- 3 x 10 22 cm 2 , i.e. compatible with 0.0. 
Note however that this effect does not influence the spectral 
index measurement. 

We have tested if the hypothesis of a constant N# agrees 
with the data : we rejected it at the 99.5 % confidence level us- 
ing a^- 2 test (^ 2 =13.41/3 d.o.f., using lcr errors in the^- 2 com- 
putation) and excluding from the fit the bin showing an increase 
of Nh at t=420 s. Because of the large redshift, we cannot set 
a meaningful constraint on the host column density with the 
PC mode restricted to energies above 0.5 keV, as already noted 
by Watson et al. (2006). Thus, in order to better constrain the 
spectral index, we fixed the extragalactic column density value 
to zero for the PC mode (after 596 seconds). Results from this 
spectral analysis are listed in Table 1 . These results are in agree- 
ment with those published by Watson et al. (2006). We note 
however that these authors reported an increase of the Nh for 
t > 9000 seconds. An inspection of the residuals indicated a 
very poor fit between 2000 and 10000 seconds (jf 2 = 2.04, 42 
d.o.f.). Adding some absorption locally to the burst fails to im- 
prove significantly the fit (xl = 2.03, 41 d.o.f.), and resulted 
only in an upper limit (Nh < 1 .2 x 10 23 crrr 2 ; however, this up- 
per limit is consistent with the value reported by Watson et al. 
2006). This poor fit could be due to some spectral component 
of the flares occurring between 2000 and 10000 seconds, thus 
we do not investigate it further. 

We restricted the spectral analysis to the early data (f <2000 
s). Figure 2 shows the variations of the softness ratio, the spec- 
tral index and the extragalactic column density as a function of 
time. The softness ratio is computed by dividing the counts de- 
tected in the 0.5-2.0 keV band by those detected in the 2.0-10.0 
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Fig. 1. Multi-wavelength light curve of GRB 050904. We present the 
X-ray data (black points) taken in the 2.0-6.0 keV band with filled cir- 
cles (WT mode data) and open circles (PC mode data), the TAROT I T 
band data with blue diamonds, the H data with solid purple squares 
(from Haislip et al. 2006; Tagliaferri et al. 2005) and the J data with 
open red symbols (from Haislip et al. 2006; Tagliaferri et al. 2005). 
For clarity, offsets by +6 and -1.5 magnitudes have been applied to 
TAROT data and H band data respectively. We indicate on the figure 
the intervals excluded for the temporal fit of the X-ray data (solid hor- 
izontal arrows). The vertical dashed line indicate the position of the 
transition between the plateau and the first flare. See electronic ver- 
sion for colors. 



keV band (i.e. the soft band is affected by absorption, while the 
hard band is not affected). We observe a global softening of the 
spectrum with time. Figure 2 also shows the variations of the 
energy spectral index, which is constant before the first flare, 
increases at the start of this flare, and then remains constant. 
A spectral fit on all the WT data before and after the start of 
the flare indicate a power law energy index of 0.23 + 0.05 and 
0.57 + 0.06 respectively (xl = 1.2, 288 d.o.f.). This latter value 
is consistent with the one derived from the following PC data 
(0.53 + 0.09). 

3.2. Optical data 

TAROT (Telescope a Action Rapide pour les Objets 
Transitoires) is a robotic observatory designed for GRB early 
detections (e.g. Klotz et al. 2005a). 

The Gamma-ray burst Coordinate Network (GCN) dis- 
tributed the SWIFT alert 81s after the beginning of the prompt 
emission of GRB 050904. The TAROT first exposure begun 
five seconds later. Twenty two unfiltered CCD images were 
taken between UTC 01:53:10.2 (t„,^+86.2s) and 02:08:09.6 
(t«r;g+ 1666.4s). Two technical problems occurred during the 
record: 

i) A bug of the scheduling software lead to stop observations 
in the range t /ng +254 to +3 12s. 
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Fig. 2. First panel (Top) : X-ray light curve of the early afterglow of 
GRB 050904 in the 2.0-6.0 keV band. We have converted the count 
rate in flux using the conversion factors derived from our spectral anal- 
ysis. Second panel) : The softness ratio of the X-ray afterglow (0.5-2.0 
keV band versus the 2.0-10.0 keV band). We corrected the bin corre- 
sponding to the flare maximum (affected by an instrumental problem) 
by substituting the observed counts in the 0.5-0.6 keV band by the 
count rate predicted by the spectral model in the same band. The cor- 
rected value is indicated by the gray symbol. Third panel) : The en- 
ergy index of the afterglow spectra. We assumed for each spectra an 
absorbed power law model. Last panel (Bottom) : the extragalactic 
column density obtained from the fit to the data. Note that the current 
calibration of the PC mode does not allow us to fit this quantity in the 
PC spectra. The symbols in panels b), c), and d) refer to the temporal 
slices listed in Table 1. 

Table 1. X-ray (0.3 or 0.5 keV to 10.0 keV band) spectral fitting re- 
sults. The xl value of the fit is 1.26 for 463 degrees of freedom. The 
Extragalactic absorption is given at z = 6.3. 



Temporal slice 


Extragalactic 


Energy index 




absorption 


a 


(s) 


(10 22 cm" 2 ) 




169-202 


18 ±5 


0.2 ± 0.2 


203-264 


4 7 +33 


0.2 + 0.1 


264-384 


< 1.93 (o.rjf ) 


0.24 ± 0.06 


384-463 


2-4 + 24 


0.65 ±0.1 


464-581 


< 1.50 (0.0 +L5 °) 


0.54 ± 0.06 


596-979 




0.53 ± 0.09 


986-1666 




0.7 ± 0.2 


2000-10000 




0.73 ± 0.08 


10000-20000 




0.74 + 0.07 


> 20000 




0.82 ± 0.08 



ii) The CCD camera cooler stopped due to a wire that broke 
during the fast slew of the telescope. As a consequence, the 
temperature of the CCD remained at -50°C until t (n ^+260s 
and then increased continuously to +14°C at t (n ^+1500s. 
Frames taken after t /ng +590s are strongly affected by ther- 
mal noise. 



We co-added frames into 6 final images in order to increase 
the signal to noise ratio. 

Table 2. Reference star characteristics nearby the position of 
GRB 050904. 



2MASS 


I 


V-R 


V-I 


Sp. Type 


013.693725 


11.68 


0.70 


1.31 


K4V 


013.732724 


12.29 


0.49 


1.01 


K2V 


013.782866 


12.39 


0.35 


0.74 


GOV 



Calibrations of TAROT data are difficult because images 
were acquired with no filter, and the high z of the GRB induces 
a non standard spectrum over the optical wavelengths. First, 
we choose three stars separated by less than 5 arcmins from the 
GRB and we calibrated their VRI magnitudes with the 80 cm 
telescope at the Observatoire de Haute-Provence, using stars 
of the SA92 field of calibration (LONEOS 2003). We com- 
pleted this study with spectra taken with the same telescope 
equipped by a low resolution spectrograph (R=150), and fit 
the star continua. Results of photometry, spectrometry and data 
from 2MASS catalog allowed us to attribute accurate magni- 
tudes and spectral types (see Table 2). We computed the so 
called superstar as the sum of the three reference stars using 
the stellar spectral flux library provided by Pickles (1998). For 
each TAROT image, we computed the Point Spread Function 
(PSF) of the superstar. We fit the optical transient (OT) of GRB 
050904 by this PSF searching the a coefficient corresponding 
to the minimum value of the (OT - ax PSF) 2 . To calibrate the 
flux of the OT, we assumed a spectral index /?=- 1.2 and a sharp 
cut-off at /l=8862 A. We modeled the response of the system 
(atmosphere, optics and CCD) and we applied it to the theoret- 
ical spectra of the OT and to the superstar. Then, we adjusted 
the flux of the OT at /l=9500 A to obtain the same a coefficient 
than that obtained from TAROT images (we call this "filter" 
the Ij band). Results are reported in Table 3, and presented to- 
gether with the X-ray data in Fig. 1. As one can see, and as 
already reported in Boer et al. (2006), TAROT has observed a 
flare in the I T band which is coincident with one X-ray flare. 

Table 3. Flux derived from TAROT images. In parenthesis, upper lim- 
its computed from PSF fit using a signal-to-noise (S/N) limit 3. 



t, rig range (s) 


Flux at 9500A (ergs cirT 2 s" 1 A" 1 ) 


S/N 


86 - 144 


(7.5xl0- 16 ) <5.3xl0" 15 


<2 


150 - 253 


4.6-7.4xl0- 15 


3.7 


312-443 


4.2-7.4xl0~ 15 


3.1 


449 - 589 


1.3 - 2.1 xl0~ 14 


4.0 


595 - 978 


(4.2xl0-' 5 )<4.2xl0- 15 


2.7 


985 - 1666 


(2.5xl0- 15 )<3.1xl0- 15 


<2 



We completed these optical observations with the ones re- 
ported in Tagliaferri et al. (2005) and Haislip et al. (2006) (see 
Fig. 1), where a steep-flat-steep evolution was found. From an 
initial decay of 1.36 +0.07 between 0.125 and 0.5 days (Haislip 



Gendre et al.: Multi wavelength observation of GRB 050904 



5 



et al. 2006), the light curve flattens to a decay of 0.7 + 0.2. At 
2.6 +1.0 days, the light curve steepens to a decay of 2.4 + 0.4, 
as usually observed in case of a so-called jet break (Rhoads 
1997). 

4. The spectral evolution of the X-ray afterglow 

The hard y-ray emission is observed up to the start of the first 
X-ray flare (t ~ 420 s, Sakamoto et al. 2005). Moreover, the X- 
ray emission before the first flare (t = 169-400 s) has temporal 
(S = 2.45+0.2) and spectral (a = 0.2 + 0.1) indexes compatible 
with the tail of the prompt one, i.e. with off-axis emission (clo- 
sure relationship 6 = 2 + a, Kumar & Panaitescu 2000), and 
the X-ray spectral index observed before the start of the flare 
is comparable to the gamma-ray one (0.2 + 0.1 vs 0.34 + 0.06 
respectively). Thus, this part of the observations may be linked 
to the prompt emission. During this phase (t<~ 400s), a signif- 
icant progressive softening is observed, correlated with a de- 
crease of the Nh value. In fact, the spectral index does not vary 
during that time (see Fig. 2). One may thus explain the soften- 
ing as the consequence of the observed decrease of the column 
density, rather than the consequence of the well known hard- 
to-soft evolution of the prompt emission. Such a decrease has 
been already observed in previous bursts, such as GRB 050730 
(Starling et al. 2005), GRB 000528 (Frontera et al. 2004) and 
GRB 980329 (Frontera et al. 2000; Lazzati & Perna 2002) al- 
though, in the latter two cases, the available statistics were too 
low to exclude a constant Nh- The decrease of the measured 
equivalent column density that we observe in the case of GRB 
050904 may be the evidence of a progressive photo-ionization 
of the initially cold gas in which the GRB occurs by the burst 
itself (Perna & Loeb 1998). 

A significant softening is observed at the start of the first 
flare; the softness ratio remains constant throughout the flare 
except at its maximum (~ 460 s) where its value changes from 
1.95 to 1.55+0.1 (lcr error). This latter variation may be due to 
a calibration problem (see Sec. 3.1 .2). Correcting for this effect 
(see Fig. 2), we obtain a softness ratio of 1.67 + 0.14 (lcr error), 
consistent with no spectral variation. 

At t=582 s, another hardening is observed. Very surpris- 
ingly, the spectral index does not vary and no excess of absorp- 
tion is detected at that time. However, the spectral fit of this part 
of the data is very poor (xl = 2-1, 13 d.o.f.). A fit with a two 
component model (e.g. a hard power law plus a soft power law, 
as expected if we observe the afterglow together with a signifi- 
cant signal from the tail of the prompt emission) is even worse 
(xl = 2.4, 1 1 d.o.f.). As can be seen from the top of Fig. 3, the 
main discrepancy arises from a bin at about 3.3 keV. We have 
tried to take this discrepancy into account by adding a narrow 
line to our spectral model rather than simply ignoring the spec- 
tral bin. We obtain a better agreement with the line parameters 
E = 3.46 + 0.15 keV and cr < 0.46 keV ( X l = 1-34, 10 d.o.f.). 
The continuum parameter changes to a = 0.66 + 0.17. Using 
this continuum model, we match a softness ratio of 1 .9, com- 
patible with the hypothesis that the softness ratio is constant 
after the first flare. The observed change of the softness ratio at 
t=582 s in Fig. 2 may thus be related only to that discrepancy 
at high energy. The statistic does not allow us to discriminate 
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Fig. 3. The X-ray spectrum of GRB 050904 integrated between 582 
and 978 seconds after the burst. Top : fit with a single power law. 
Bottom : fit with a power law continuum and a Gaussian line. 

between a spurious deviation and a real emission line. This de- 
viation could be due to instrumental problems (such as back- 
ground fluctuations) rather than a real physical process. 

After the initial 2 ks of data, the XRT light curve features 
several flares (see Fig. 1). X-ray spectra extracted during the 
3rd, 4th and 5th flares indicate a small softening of the spectral 
properties. We cannot investigate the underlying decay law due 
to the large flaring activity. 

5. The stellar progenitor and the surrounding 
medium of GRB 050904 : observation of a 
termination shock 

We have used the closure relationships given by Chevalier & Li 
(2000), Sari et al. (1998), Sari & Piran (1999) to investigate 
the surrounding medium of GRB 050904 and on the fireball ge- 
ometry. We considered the X-ray data between the end of the 
first flare (t~580 s) and ~ 1666 s. In the optical band, we sepa- 
rate the data between early and late ones, as observed between 
0.1 and 0.5 days, and between 0.5 and 2.6 days respectively. 
The transition date (0.5 days) is chosen because both Haislip et 
al. (2006) and Tagliaferri et al. (2005) observe a flattening at 
that time. The early optical data occur during the 4th and 5th X- 
ray flares. We thus prefer to conservatively discard these data, 
and use only the late ones for the closure relationship study, 
adopting the a and 6 values given by Tagliaferri et al. (2005). 
We list the closure relationship results in Table 4. 

The X-ray data are marginally consistent with the jet hy- 
pothesis. However, in this case, one should not expect any fur- 
ther steepening in the light curve (Rhoads 1997), contrary to 
that observed by Tagliaferri et al. (2005) at about 2.6 days. 

An ISM scenario (slow cooling phase with v e above the 
X-ray band) is only in marginal agreement with the early X- 
ray data (t<2000s, see Table 4). Moreover, the expected differ- 
ence between the X-ray and optical light curve decay indexes 
should be 6 X - 6 — or 5 X - 6 — 0.25. The observed value is 
5x - S = 0.8 + 0.4, again only marginally consistent with the 
expectations. On the other hand, the late broad-band (radio-to- 
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Table 4. Closure relationships in the standard fireball model computed using the spectral and temporal information. We give the relationships 
for the J band after 0.5 days. The specific frequency is v c (the cooling frequency) in the slow cooling regime, or v m (the injection frequency) in 
the case of fast cooling. All errors are quoted at the 90 % confidence level. 



Medium class Cooling Specific Closure Expected X-ray late 

and regime frequency relationship value optical 



geometry 




position 








(0.5-10.0 keV) 
(582-1666 s) 


(7 band) 
(0.5-2.6 days) 


Isotropic Wind 


Fast 


v m < V 


6- 


1.5a 


-0.5 


0.6 ± 0.3 


-1.1 ±0.6 






Vm > V 


6- 


0.5a 


0.0 


1.2 ±0.3 


0.1 ±0.4 




Slow 


V c < V 


6- 


1.5a 


-0.5 


0.6 ± 0.3 


-1.1 ±0.6 






V c > V 


6- 


1.5a 


0.5 


0.6 ± 0.3 


-1.1 ±0.6 


Isotropic ISM 


Fast 


Vm < V 


6- 


1.5a 


-0.5 


0.6 ± 0.3 


-1.1 ±0.6 






Vm > V 


6- 


0.5a 


0.0 


1.2 ±0.3 


0.1 ±0.4 




Slow 


V c < V 


6- 


1.5a 


-0.5 


0.6 ± 0.3 


-1.1 ±0.6 






v c >v 


6- 


1.5a 


0.0 


0.6 ± 0.3 


-1.1 ±0.6 


Jetted fireball 


Slow 


v c <v 


5- 


-2a 


0.0 


0.4 ± 0.3 


-1.7 ±0.8 






V c > V 


5- 


-2a 


1.0 


0.4 ± 0.3 


-1.7 ±0.8 



X-rays) observations clearly agree with an ISM environment 
having a very high density and v c below the optical band, as 
recently shown by Frail et al. (2006). This model predicts that 
the cooling frequency is below the optical band even at early 
(~ 500) times. However, its extrapolation in the X-ray band 
at these times (black-dotted line in Fig. 4) cannot describe the 
data: the predicted flux is at least one order of magnitude be- 
low the observed one and the extrapolated light curve is flatter 
than the observed best fit decay (red-dash-dotted line), as al- 
ready anticipated above. Using the closure relationships, the 
hypothesis that the interstellar medium explains the X-ray data 
is excluded at the 5.8<x level. We thus conclude that while an 
ISM describes well the late broad-band observations (Frail et 
al. 2006), it does not reproduce satisfactorily the early X-ray 
data (t<2000s). 

A wind environment is clearly favored by the early 
(t<2000s) X-ray observations (slow cooling, v c above the X- 
ray band, see Table 4). In a wind medium v c increases with 
time (Chevalier & Li 2000), thus we expect also the optical 
band to be below this frequency at late times. However, this 
hypothesis is rejected at the 4.2<x level. A similar case with a 
flat spectral index and a typical temporal decay was observed 
only for GRB 040106, and attributed to a wind density profile 
(Gendre et al. 2005). Combined with the results by Frail et 
al. (2006), this suggests that the fireball could be expanding 
into a medium whose density profile at small distances is that 
expected for a wind (« oc r~ 2 ), and becomes constant at larger 
radii. 

In this hypothesis, a wind termination shock is expected 
to mark the transition between the two different environments. 
The termination shock radius, R t , is defined as a function of the 
wind density parameter A, and the ISM density n (Chevalier et 
al. 2004; Panaitescu & Kumar 2004): 

R, = 1.1 x 10 18 A°V 05 (1) 

where R t is in units of cm. The termination shock crossing 
time can be directly estimated from the observations to be lo- 
cated between 0.019 days (the latest time at which the X-ray 
closure relationships can be computed and favor a wind envi- 



ronment) and 0.5 days (the earliest time at which the optical 
data favor an ISM). One can then use the standard fireball pre- 
scriptions that give the fireball radius as a function of time, 
initial energy £52 (in the remaining we will use the standard 
notation A x = 10 A A) and A* to derive an independent estimate 
of the termination shock radius as a function of these two quan- 
tities (Chevalier et al. 2004): 

^.ssxio'^fA; - 5 ^) ' 5 (2) 

where t,j is the shock crossing time in units of days. 
Combining equations (1) and (2), one can derive the possible 
values of A„: 

£° 2 5 n a5 7.24 x 10~ 3 < A* < £° 2 5 « 05 3.69 x 10~ 2 (3) 

Using the best fit results by Frail et al. (2006), n=680 and 
£52=88, we obtain 1.77 < A„ < 9.03. 

We thus derive the position of the termination shock radius, 
R„ to be 0.018 pc < R, < 0.041 pc. Note that this burst is very 
energetic (Kann et al. 2006) and surrounded by a very dense 
medium (Frail et al. 2006), and thus the termination shock po- 
sition of GRB 050904 may not be typical. However, this value 
is compatible with the upper limits obtained by Chevalier et al. 
(2004) on several bursts. 

Using the prescriptions given by Yost et al. (2003), we 
have checked if there is a solution compatible with the X- 
ray data within the range of values for A, derived using the 
best fit by Frail et al. (2006). As shown in Fig. 4, setting 
A, = 1.8, we nicely describe the early X-ray data (black solid 
line before 2000 s). This implies a termination shock radius 
of R, ~ 1.8 x 10~ 2 pc. Note that we find e B = 10~ 6 and e e = 
8.5 xl0~ 2 , that are different from Frail's values. However, such 
a change is reasonable because of the shock presence. In fact, 
the microphysic parameters are allowed to vary during the fire- 
ball evolution even if there is no termination shock (Panaitescu 
et al. 2006). Thus, a wind plus ISM scenario can provide a 
global model for the broad-band observations. 

With these parameters, we find a wind density of ~ 170 
cm" 3 at the termination shock radius, giving a density jump of 
~ 4 as implied by eq. 1. 
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100 1000 10 4 10 5 10 6 

Time (s) 

Fig. 4. Global modeling of the broad-band data using a wind (A, = 1.8, e B = 10~ 6 , e e = 8.5 x 10~ 2 ) plus ISM (n = 680, e s = 0.015, e e = 10~ 2 ) 
scenario (solid lines). We present in black the X-ray data and in green the optical data (circles are TAROT data rescaled in the J band, squares 
are data taken from Tagliaferri et al. 2005; Haislip et al. 2006, rescaled by a factor 0.01 for clarity purpose). The red-dot-dashed line is the 
best fit decay law to the X-ray data between ~ 600 s and ~ 1700 s. The dotted lines represent the extrapolation of the ISM model by Frail et 
al. (2006) to early times. The dashed lines are the extrapolation of the wind model to late times. We assumed an instantaneous wind-to-ISM 
transition at t ~ 1700 s (see text), smoothed by the effect of the fireball curvature (At ~ 3400 s). See electronic version for colors. 



Finally, an order of magnitude estimate of the intrinsic ab- 
sorption implied by the wind environment is given by Piro et 
al. (2005): 

N H ,22 = 3— (4) 
13 

where ro is the radius at which the prompt emission is pro- 
duced. Using A» = 1.7 and assuming 03 = 1, we find con- 
sistency with the value we observed at the start of the X-ray 
observation. 

6. Multi-wavelength modeling of the first flare 

With the Swift observations, it appears clear that flares are very 
common in GRB light curves during a time interval that goes 
from hundred to thousand of seconds (O'Brien et al. 2006). 
The evidence of such a common flaring activity has favored 
the development of several models that attempt to explain the 
nature of X-ray flares, both in the framework of Internal Shock 



and External Shock. In the following, we comment on the first 
flare observed simultaneously in the X-ray and optical bands in 
light of these models. 

6.1. Delayed external shock scenario 

The X-ray spectrum of the first flare is softer than the preced- 
ing plateau and is consistent with the following afterglow emis- 
sion. Similar properties have been previously observed in other 
bursts (e.gGRB 01 1 121, XRR 01 121 1 and XRF 01 1030, Piro 
et al. 2005; Galli & Piro 2006), and have been explained in the 
framework of a delayed external shock, expected e.g. in a thick 
shell fireball scenario (Sari & Piran 1999). In this framework 
both the X-ray and the optical flares are produced by the ex- 
ternal shock, thus explaining the coincidence between the two 
flares. 

Very recently, Lazzati et al. (2005) have shown that a cen- 
tral engine releasing most of its energy reservoir during the fi- 
nal stages of its activity (t eng ), implies that the afterglow emis- 
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Fig. 5. X-ray light curve of GRB 050904 produced by a thick shell 
fireball expanding in a wind with the onset of the afterglow shifted 
to to = 464 s. The model parameters are E 52 =200, r =200, A„=0.5, 
e e =0.008, e B = \Q-\p=2A. 

sion is described by a simple power law only if the origin of 
the time to is of the order of t eng . This implies that the begin- 
ning of the afterglow emission to is the time of flare appearance 
t flan (Galli & Piro 2006). We used the model of Galli & Piro 
(2006), that take this effect into account, to check if it can ex- 
plain the multi-wavelength flare. 

This model is a single component model that does not take 
into account the contribution of the prompt emission (that could 
affect the rise of the flare), nor the contribution from the reverse 
shock (that does not cross instantaneously the thick shell). 
Thus, it can fail to describe the rising of the flare emission, 
but it can be applied to the decay of the first flare and the fol- 
lowing X-ray data (excluding all subsequent flares). In the case 
of a fireball expanding in a wind with v < v c , a — (p - l)/2 
(Chevalier et al. 2004). We derive a ~ 0.55, which implies 
p = 2.1. We have fixed £52=200. The remaining model pa- 
rameters, i.e. A*, € e and eg are determined in order to describe 
the data. We find a solution with parameters Fq=200 (where To 
represent the initial Lorentz factor), e e =0.008, e B = 10~ 6 and 
A, = 0.5 (see Fig. 5). This value of A* implies 1 < n < 20, 
compatible with the range 10 < n < 10 4 given by Frail et al. 
(2006). 

However, some inconsistencies are implied by this sce- 
nario. First of all, in a thick shell, the afterglow peak is at the 
time when the prompt emission ceases (Sari & Piran 1999). 
As discussed in Sect. 4, the data preceding the flare are consis- 
tent with high latitude emission observed after the end of the 
prompt phase. Under this assumption, there is a temporal gap 
between the end of the prompt emission and the flare. Thus, 
such a gap is not compatible with relating the appearance of 
the flare with the onset of the afterglow from a thick shell. 

Moreover, even without considering the data before the 
flare, no solution can reproduce the optical observations in a 
wind case : the calculated optical flux at the time of the flare 
is ~ 10~ 26 erg ctrT 2 s _1 Hz T 1 , about 1.5 order of magnitude 
below the measured flux. In fact, with the cooling frequency 
v c above the X-ray band, the predicted broad band optical to 



X-ray spectral index value (a ox ) is ~ 0.6. Instead, connecting 
the mean flux at 5 keV and the optical flux at the time of the 
flare, a ox is ~ 1 . 1 . As a consequence the extrapolation of the 
X-ray flux in the It band falls below the optical data. A value 
of a ox ~ 1.1 at the time of the flare can be obtained assuming 
that v c is below the optical range during the flare, with p — 2.1, 
but this is not compatible with the closure relationships. 

6.2. Possible two-component scenarios 
6.2.1 . The role of the reverse shock 

GRB 050904 is one of the few GRBs (the others are GRB 
990123, GRB 041219A, GRB 050401, GRB 06011 IB, GRB 
060124, GRB 060904B) for which a prompt optical emis- 
sion was observed simultaneously with the high energy one 
(Akerlof et al. 1999; Fan & Wei 2005; Rykoff et al. 2005; 
Boer et al. 2006; Klotz et al. 2006a; Romano et al. 2006; 
Klotz et al. 2006b; de Ugarte Postigo et al. 2006). Boer et al. 
(2006) underlined the similarities between GRB 990123 and 
GRB 050904 in terms of optical flare brightness (once account- 
ing for the difference in the distance), and broad-band spectrum 
of the prompt emission : in both cases the extrapolation of the 
high energy data to the optical band falls well below the optical 
points (Corsi et al. 2005; Boer et al. 2006), suggesting a differ- 
ent origin for the low and high energy emission mechanisms. 
In the case of GRB 990123, this hypothesis was further sup- 
ported by the lack of coincidence between the peak observed 
in the optical and those observed in the high energy light curve 
(Corsi etal. 2005). 

In the standard fireball model (Rees & Meszaros 1992; 
Meszaros & Rees 1997; Panaitescu et al. 1998), when a rel- 
ativistic ejecta moves into the cold ISM, two shocks form, an 
outgoing one that propagates into the ISM (the Forward Shock, 
FS) and a Reverse Shock (RS) that propagates into the ejecta 
(Sari & Piran 1999). 

The typical synchrotron frequency of the RS emission is 
in the IR-to-optical region. If the signal observed in the low- 
energy part of the prompt spectrum is related to the RS, while 
prompt y- and X-ray emission to internal shocks (or late inter- 
nal shocks for the late time X-ray flares), one should not expect 
the low energy emission to fall on the extrapolation of the high 
energy one; moreover, there should not be any correlation be- 
tween the pulses observed in the high energy light curve and the 
peak time of the optical flare. This is in fact what was observed 
for GRB 990123. In the case of GRB 050904, the high energy 
spectrum cannot be extrapolated at low energies, but the optical 
flare observed by TAROT appears to peak simultaneously with 
the first X-ray flare. In the internal shock-RS scenario, such a 
coincidence should be fortuitous. 

Recently Fan & Wei (2005) explored the possibility of re- 
lating X-ray flares to RS synchrotron emission appearing above 
the FS one, when the physical parameters (e e and eg) in the RS 
and the FS are different. In this framework, one could think of 
relating both the optical and the first X-ray flare to RS emission. 
The temporal coincidence between those two flares would be a 
natural expectation of this model. However, Wei et al. (2006) 
noted that even in this scenario, the X-ray spectrum should be 
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a power-law extension of the optical emission, but the obser- 
vations show that the optical-to-X-ray spectrum cannot be de- 
scribed by a simple synchrotron spectrum, even taking into ac- 
count the possible presence of a spectral break between the op- 
tical and the X-ray band. 

6.2.2. Flares in a RS synchrotron plus inverse 
Compton scenario 

Recently, Kobayashi et al. (2005) suggested the possibility of 
explaining X-ray flares via synchrotron self-Inverse Compton 
(IC) radiation from the RS. The prompt optical flare should be 
associated with synchrotron emission from the RS while the 
simultaneous X-ray flare should be the result of synchrotron 
photons being up-scattered via the IC process in the X-rays. 
The RS evolution depends on the initial Lorentz factor of the 
shell itself (Sari 1997). In the so-called thin shell case, the RS 
remains Newtonian during all the crossing time and the effec- 
tive energy extraction takes place at the deceleration time td- 
On the other hand, if the initial Lorentz factor of the shell is 
sufficiently high, i.e. greater than the critical value F c : 



T c ~ 130 



'1 +z 



3/8 



p l/8 T -3/8 -1/8 
n 52 1 2 " 



(5) 



then the RS becomes relativistic while crossing the shell 
(Kobayashi et al. 2005). In Eq. 5, Ti is the burst duration. 
This corresponds to the thick shell case. In a thin shell fireball 
expanding in an ISM, both the low and high energy emissions 
should thus peak at tj and the temporal coincidence between 
the optical and the X-ray flares would be a natural expecta- 
tion for this model too. In particular, in the thin shell case, the 
peak time of the X-ray and optical flares (~ 460 s) constrains 
the model parameters according to the following relationship 
(Kobayashi et al. 2005): 



(?) 



•1/3 



460 s 



(6) 



Setting that £52 ~ 200, z ~ 6.3, for n ~ 100 it is T ~ 120 
and F e ~ 170 for T ~ 225 s, consistent with the hypothesis of 
being in a thin shell case (i.e. To < F c ). 

Using an additional emission mechanism like IC gives a 
natural explanation for the impossibility of extrapolating the 
spectrum of the X-ray flare down to the optical data points. 
In particular, we know that the 0.5-10 keV spectrum of the 
flare is rather flat, with a spectral index of a ~ 0.6. Setting 
p ~ 2.1, the spectrum of the flare can be explained assuming 
that the X-ray band is between v^f = 2y 2 v m and v ! c c = 2y L c v c 
around the deceleration time (Sari & Esin 2001). Considering 
that y m ~ 100, if v^ c ~ 1.0 keV and v'f ~ 10 keV, then v m ~ 
10 13 Hz and v e ~ 3 x 10 13 Hz, that are reasonable values for 
the RS synchrotron break frequencies at the deceleration time 
(Kobayashi et al. 2005). With p = 2.1, the RS synchrotron 
flux at ~ 1 keV is ~ 1/1060 times the optical one. The ratio 
between the peak flux observed by TAROT between 449 s and 
589 s, and the mean flux observed in the same temporal bin 
at 1 keV is ~ 1/750 (i.e. a factor of ~ 1.4 greater than that 
extrapolated from the RS synchrotron emission), while the ratio 



between the same optical flux and the peak flux of the X-ray 
flare at 1 keV is ~ 1 /280 (i.e. a factor of ~ 4 greater than 
that extrapolated from the RS synchrotron emission). Thus, to 
explain the broad-band spectrum, the RS IC peak flux should be 
a factor of ~ 1 .4 - 4 greater than the RS synchrotron flux at the 
same frequency (which is also compatible with the observed 
2.0 - 6.0 keV flux increase by a factor of ~ 2 in the rising 
part of the flare). As shown by Kobayashi et al. (2005), an IC 
bump as high as 6 times the synchrotron RS emission can be 
explained with a reasonable choice of parameters. 

The model proposed by Kobayashi et al. (2005) thus ap- 
pears to be a viable one to explain the broad-band spectrum 
and the temporal coincidence between the optical and the X-ray 
flares. However, one problem arises : explaining the first X-ray 
flare via IC emission with v^f < vx < v' c c , implies that the de- 
caying part of the flare should be as (f/fj)~ (3/,+1)/3 (Kobayashi 
et al. 2005), which for p = 2.1 gives 6 = 2.4, too shallow to 
agree with the observations (the X-ray emissions declines as 
r 12 6 between 464 s and 534 s). 

Moreover, in order for the RS emission to be dominant on 
the FS one, the peak frequency of the FS at the deceleration 
time should be << v^; combining this with the fact that 
the spectral index of the early X-ray afterglow is rather flat, 
it should be min(v m ,v c ) < vx < max(v m ,v e ) until ~ 10 s s. 
However, as shown in Table 4, in this regime the expected tem- 
poral decay is too flat to be compatible with the observations. 

Finally, the presence of a plateau before the optical flare ob- 
served by TAROT should also be considered: Wei et al. (2006) 
shown that a plateau plus an optical flare could be produced if 
the outflow had a more complex structure (e.g. two components 
with different Lorentz factors, widths and isotropic energies). 
However, a fit to the optical light curve gives a value for eg 
too high for being compatible with the idea of explaining the 
X-ray flare as IC emission from the RS. In fact, a high magne- 
tization suppresses the importance of IC emission with respect 
to synchrotron one. 



6.2.3. Internal shock scenario 

Wei et al. (2006) have recently shown how the problem of the 
steep temporal decay of the X-ray flare can be solved invoking 
a late internal shock model, where the optical flash comes from 
late internal shock synchrotron emission while the first X-ray 
flare is produced by late internal shock IC emission. In this 
model the temporal coincidence between the optical and the X- 
ray flare is a natural expectation as it is in the model proposed 
by Kobayashi et al. (2005) and there is no optical to X-ray 
extrapolation problem. 

Finally, Zou et al. (2005) interpreted the early-to-late 
time multi-band observations assuming that all the highly vari- 
able X-ray emission of GRB 050904 originated from internal 
shocks developing when faster shells, continuously ejected by 
the central engine, overtake a slower one. In this model, the 
main contribution to X-ray emission comes from relativistic 
RSs while the optical emission is ascribed to Newtonian FSs. 
This model is very different from the standard assumptions of 
the canonical fireball model, but seems to be able to fit the pro- 
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longed flaring activity observed in the X-rays. However, we 
notice that the optical flare and the overall decay behavior of 
the X-ray light curve still need to be modeled in detail. We also 
underline that this model assumes a super-long central engine 
activity, i.e. that the shell ejection process lasts 10 4 - 10 5 s. 



7. Conclusion 

We have analyzed multi-wavelength observations of GRB 
050904 obtained with TAROT and SWIFT. Before the first X- 
ray flare, the X-ray decay and spectral indexes are consistent 
with the hypothesis of observing the tail of the prompt emis- 
sion by curvature effect. After the flare, the X-ray emission 
is consistent with the hypothesis that the fireball expands in 
a wind medium. This would privilege a stellar progenitor. We 
estimate the wind density parameter to be 1.77 < A* < 9.03. 
After 0.5 days, the data are consistent with the hypothesis that 
the fireball is expanding in an interstellar medium. The hypoth- 
esis that a single kind of surrounding medium can explain all 
data is clearly rejected : the probability that the wind medium 
(which explains the X-ray data) can explain the optical data is 
1.7xl0~ 5 (exclusion at 4.2cr) and the probability that the in- 
terstellar medium (which explains the optical data) can explain 
the X-ray data is 4xl0~ 9 (exclusion at 5.8<x). This implies that 
the fireball may have crossed a termination shock. From the 
temporal constraints on the crossing time, we estimate the ter- 
mination shock position to be 1.8 x 10 2 pc < R t < 4.1 x 10~ 2 
pc. We observe a significant excess of absorption at the start of 
the X-ray observation, Nh = (18 ± 5) x 10 22 crrT 2 , compati- 
ble with the wind expectations. We further observed a decrease 
of this column density during the first seconds of the observa- 
tion, significant at the ~ 3<x level, that could be explained by 
photo-ionization of the medium by the burst. We investigated 
the simultaneity of the X-ray and optical flares in the context 
of several models. None of them can explain all the available 
data. In particular, the optical-to-X-ray spectrum cannot be de- 
scribed in a simple synchrotron scenario. Adding the contri- 
bution of IC emission in the X-ray band, a RS model fails to 
explain the steepness of the flare decay, while a late internal 
shock is favored. 

Due to the time dilation associated with its large distance, 
GRB 050904 offered for the first time the possibility of well 
sampling the prompt-to-afterglow transition, and the very early 
afterglow phase at several wavelengths. Its observation opened 
more questions than confirming canonical models. SWIFT is 
well suited to detect high redshift events, and hopefully other 
observations will be available to close those questions. 
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